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5-H ' Abstract The changes that vacancies produce in the properties of hep solid ^He 

are studied by means of quantum Monte Carlo methods. Our results show that 
the introduction of vacancies produces significant changes in the behavior of solid 
^He, even when the vacancy concentration is very small. We show that there is 
an onset temperature where the properties of incommensurate "^He change signif- 
icantly. Below this temperature, we observe the emergence of off-diagonal long 
range order and a complete spatial delocalization of the vacancies. This temper- 
ature is quite close to the temperature where non-classical rotational inertia has 

^ I been experimentally observed. Finally, we report results on the influence of va- 

Q . cancies in the elastic properties of hep ^He at zero temperature. 

O ■ PACS numbers: 67.80.-s,02.70.Ss,67.40.-w 



<*' • 1 Introduction 



The counterintuitive concept of supersolidity can only be understood in terms of 
extremely quantum matter. The simultaneous existence of spatial solid order, char- 
acteristic of the solid state, and superfluidity, property that in principle requires of 
^_^ ' mass movement without friction, is hardly understandable in usual solids. The 

natural candidate for this fascinating possibility has been along the time always 
the same, solid '*He. The solid phase of ^He is not a normal solid, in the classical 
meaning of this term. ''He atoms have so small mass and shallow interaction that 
its stable condensed phase, in the limit of zero temperature, is a liquid and thus a 
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finite pressure is required to crystallize it. Even in the crystal phase, the zero-point 
motion of the ^He atoms is not fully depressed as it is explicitly quantified by its 
large Lindemann ratio and excidingly large kinetic energy, with respect to typical 
classical values. Furthermore, the exchange frequency, which is absolutely absent 
in classical solids, is low but not zero. Altogether makes the study of a quantum 
solid a very interesting topic, the possibility of a supersolid scenario being one of 
its more stimulating possibilities. 

The recent experimental findings by Kim and Chan^^ on the existence of a 
non-zero fraction of non-classical moment of inertia (NCRI) in torsional oscillator 
measurements have revived this topic that emerged in the past as an hypothetical 
theoretical conjecture. This unexpected result has been corroborated by other lab- 
oratories, although the dispersion in the size of the effect is large ^. From the very 
beginning, it was clear that the way in which the crystal is produced, the annealing 
during its growth and the purity of the '^He sample were, among others, relevant 
parameters which introduce changes of the NCRI fraction of an order of magni- 
tude. The decoupling of a part of the mass, measured in the torsional oscillator, is 
observed at an onset temperature which also depends on the particular conditions 
of the experiment, but its fluctuation is sizeably smaller than the one obtained for 
the NCRI fraction (superfluid fraction). Very close to the onset temperature for 
supersolidity it has been observed an increase in the shear modulus of hep '^He, 
its temperature dependence being similar in shape to the one of NCRI^. The in- 
terplay between this anomalous elastic behavior and the NCRI effect observed in 
torsional oscillators has been object of theoretical debate: it seems clear that part 
of the NCRI effect can be attributed to elasticity but probably not all. 

To date, there is not a complete theoretical understanding of the observed phe- 
nomena. A point in which the community working in this field has reached an 
overall agreement is the absence of any supersolid signal in a perfect crystal, i.e., 
a crystal where the number of particles and the number of sites matches exactly. 
The supersolidity is then attributed to alterations of the perfect crystal due to the 
disorder introduced by defects, that inevitably appear during its growing process 
from the liquid phase. The major role has been assigned to dislocations and their 
mobility in relation with the pining/depining of -^He impurities to them. But, even 
if dislocations can be in the origin of experimental NCRI they can not be the only 
explanation since the superfluidity of hep ^He in bulk or in vycor is nearly the 
same, whereas the dislocation density is reasonably expected to be very differ- 
ent. Another possibility, that has been theoretically explored from long time is the 
presence of a fraction of point defects, i.e., vacancies. 

Vacancies were originally proposed by Andreev and Lifshitz^ as a mecha- 
nism for creating a supersolid phase in solid helium. This scenario was rebuked 
by experimental and theoretical findings showing that vacancies are energetically 
too expensive to be created by thermal activation. However, Rossi et al. ^ showed 
using a variational estimate that wavefunctions that are able to describe well the 
equation of state of solid "^He support a finite vacancy concentration between 10~ 
and 10~^. This possibility is yet to be tested with first principles calculations, as 
the number of helium atoms that corresponds to such low concentrations is still 
too large to handle. Even such low vacancy concentrations should be experimen- 
tally significant, given the observed influence of sub-ppm concentrations of ^He^-. 



Notice that as a quasiparticle in hep '*He, a vaeaney is in many ways similar to a 
^He impurity. 

Even if thermal activation of vacancies seems improbable at very low temper- 
ature it is reasonable to consider that a significant number of vacancies can be 
introduced into solid samples during their growth. Eventually, some of them can 
disappear by migrating to dislocations or grain boundaries but one cannot exclude 
a priori that a tiny fraction of vacancies can dissolve into the bulk crystal. It was, 
however, argued ^ that even if initally vacancy-rich, solid hep ^He would phase 
separate into a vacancy-rich phase and a perfect, insulating crystal and therefore 
any growth-introduced vacancies will be effectively removed from the experimen- 
tal samples. On the other hand, our detailed studies with several vacancies show 
no sign of vacancy clustering^''". 

In this work, we report recent results on the properties of vacancies in a fully 
quantum crystal like '^He obtained using different quantum Monte Carlo methods. 
In Sec. 2, we describe the microscopic methods used in the present analysis. In 
Sec. 3, we present PIMC results on the one-body density matrix of hep solid "^He 
as a function of temperature, showing the onset temperature where both Bose- 
Einstein condensation and vacancy delocalization appear. Sec. 4 comprises results 
of the shear modulus at zero temperature and as a function of the pressure. Finally, 
an account of the main conclusions of the present work is included in Sec. 5. 



2 Quantum Monte Carlo methods in the study of soUd He 

A quantitatively accurate study of solid "^He is demanding due to the high density 
of the system and strong interparticle correlations. If a microscopic approximation 
to the system is pursued, the most powerful tool is quantum Monte Carlo with sev- 
eral methods to deal properly with zero or finite temperature simulations. The final 
goal is to get relevant microscopic information, on both energy and structure, start- 
ing directly from the Hamiltonian of the system. The helium interatomic potential 
is accurately known and there are several models that are able to describe its equa- 
tion of state very well. In the present work, we have used an Aziz potential '' that 
has proven to be very accurate in the reproduction of the equation of state P{p) in 
both liquid and solid ^He. 

Calculations at zero temperature have been performed with diffusion Monte 
Carlo (DMC). DMC is a zero-temperature first-principles method which can ac- 
cess exactly the ground state of bosonic systems. It is a form of Green's Func- 
tion Monte Carlo which samples the projection of the ground state from the initial 
configuration with the operator exp [— ( Jf — Eo)t]. Here, Jif is the system Hamil- 
tonian, Eq is a norm-preserving adjustable constant and T is the variable which 
corresponds to imaginary time. The simulation is performed by advancing in T 
via a combination of diffusion, drift and branching steps on walkers (sets of 3A'^ 
coordinates) representing the wavefunction of the system. 

One of the advantages of DMC lies in the convenient incorporation of the im- 
portance sampling. The imaginary time evolution of the walkers is "guided" dur- 
ing the drift stage by a guiding wavefunction ^q, which is usually a good guess 
for the wavefunction of the system. When a guiding wavefunction is used, the ex- 
pectation value of an operator jz/ computed with DMC results in a value equal to 
(00 \£/\ (J)g), where 00 is the ground state wavefunction of the system. This leads 



to a common misunderstanding regarding the bias resulting from the choice of 
^G- In fact, it is straightforward to show that for the Hamikonian Jf and any 
operator commuting with it, the expectation value is computed exactly within sta- 
tistical error. That is, if [^, Jf ] = 0, then {^q \j2/\ ^q) = (0o |=e/| (/)o) for any ^q 
such that ((/)o|(/>g) 7^ 0. In practice, one usually demands that the DMC results are 
not sensitive to small changes in parameters describing (J)g, if the parameters are 
sufficiently close to the optimal values. The use of the importance sampling is the 
reason why DMC results agree so precisely with a wide range of energetic and 
structural experimental data, both for liquid and solid ^He. 

The wavefunction that we use for the importance sampling of solid helium 
is a symmetrization of the well-known Nosanow-Jastrow— wavefunction. The 
symmetrized version has the form 



/iVp \ /W. Wp \ 

0sN.=(n/(i'-'-'-./i)j (ni^(i'-<-'^i)j- (1) 



The first term in (/)snj is the McMillan product of pair correlation functions f{r) = 
exp [—l/2{b/ry] . The second term describes the lattice structure. Here r and / 
describe respectively positions of the A^p atoms and A^s lattice sites. The localizing 
function g{r) is set to a Gaussian, g{r) = exp( — 1/2//-^). The two parameters, b 
and y, are obtained from optimization of the unsymmetrized form of the function, 
as described in Ref. '^. This allows for excellent convergence, including insensi- 
tivity of results with respect to changes in the parameters, as described above. 

The symmetrized wavefunction of Eq. ([TJ is an excellent and, at the same 
time, simple description of a quantum crystal. Already on the variational level, 
this wavefunction is capable of predicting the solid-to-liquid transition for ^He 
(optimized 7 experiences a discontinuous jump to y — 0, corresponding to the 
liquid phase), which by itself is a very remarkable result. The second lattice term 
in (/)sNj is a product of sums over all particles. The product is maximized when 
the smallest possible number of the sums is vanishingly small. This corresponds 
to largest possible number of lattice sites being occupied. The product therefore 
includes all possible permutations of atoms on the lattice, but the atoms are more 
likely to occupy maximum number of sites in each permutation. Interstitials are 
allowed but have a built-in probability penalty, as it should be. Exchanges between 
atoms are allowed by construction. Vacancies may be introduced by using A^p < 
A's, which only changes the norm of (/>sni- Structural defects may be created by 
perturbing the set of the lattice sites {I}. This wavefunction was first introduced 
and tested with DMC for quantum solids by Cazorla et al.— and has been since 
then applied to a range of problems including elasticity— and point defectsi^^ in 
solid ^He, and even quantum crystals of atoms excited to the Rydberg states^. 

We are also interested in the behavior of vacancies in solid ^He at finite tem- 
perature. In this case, the appropriate theoretical tool is the path integral Monte 
Carlo (PIMC) method. Starting from the Hamiltonian J^ and the temperature 
T — (ksP)'^ of the system, it is possible to rewrite the partition function Z mak- 
ing use of the convolution property that the thermal density matrix p (R' , R; /3 ) = 



(R'|e-'^^|R) satisfies, 



M 



/iVl 
ndRiP(Ri,Ri+i;e) , (2) 

i=i 

with e = P/M, and the boundary condition Rm+i = R\- 

PIMC describes the quantum A^-body system considering M different config- 
urations Ry of the same system, whose sequence constitutes a path in imaginary 
time. This means that the A'^-body quantum system is mapped onto a classical sys- 
tem of A' ring polymers, each one composed by M beads. The different beads can 
be thought as a way to describe the delocalization of the quantum particle due to 
its zero-point motion. 

For sufficiently large M, we recover the high-temperature limit for the thermal 
density matrix, where it is legitimate to separate the kinetic contribution from the 
potential one (primitive action). In this way, it is possible to reduce the systematic 
error due to the analytical approximation for p below the statistical uncertainties 
and therefore to recover "exactly" the thermal equilibrium properties of the sys- 
tem. However, the primitive action is too simple for studying extreme quantum 
matter and a better choice for the action is fundamental to reduce the complexity 
of the calculation and ergodicity issues. Using the Chin action '^■'^, we are able to 
obtain an accurate estimation of the relevant physical quantities with reasonable 
numeric effort even in the low temperature regime, where the simulation becomes 
harder due to the large zero-point motion of particles. 

An additional problem we have to deal with when approaching the low tem- 
perature limit with PIMC simulations arises from the indistinguishable nature of 
'^He atoms. In the path integral formalism, the exchanges between L different par- 
ticles are represented by long ring polymer composed by L x M beads. If we study 
a bosonic system, the indistinguishability of the particles does not affect the posi- 
tivity of the integrand function in Eq. |2]and thus the symmetry of Z can be recov- 
ered via the direct sampling of permutations between the ring polymers. A very 
efficient sampling scheme that we have used in the present study is provided by 
the Worm Algorithm '^. 

The formation of long permutation cycles is a frequent event at very low 
temperature and polymers which close themselves winding the periodic bound- 
ary conditions of the simulation box can appear in the PIMC configurations. The 
winding number, that is the net number of times the paths of the polymers wind 
around the periodic cell, is an important quantity since it is related to the superfluid 
fraction of the system^". In particular, the appearance of polymers presenting non- 
zero winding numbers in PIMC configurations gives indication of superfluidity in 
the simulated system. 



3 Vacancies at finite temperature 

The presence of vacancies in solid '^He produces at zero temperature the emer- 
gence of supersolidity, even when the vacancy concentration is very small. The 
influence of vacancies in solid ^He is reasonably expected to change with temper- 
ature, but an accurate estimation of temperatures at which supersolidity appears 
was lacking until recently. Our aim has been to gain a deeper understanding of this 



behavior by performing comprehensive PIMC simulations. It is worth mentioning 
that the efficiency of the sampling in PIMC, when the temperature approaches the 
zero limit, drops progressively by the low acceptance ratio of the sampling move- 
ments. In order to overcome, at least in part, these technical issues it is crucial to 
work with an accurate action, that allows for reducing the number of terms (beads), 
and a good permutation sampling as the one provided by the worm algorithm. 

We have studied the properties of incommensurate solid ^He at finite temper- 
ature carrying out PIMC simulations of N = 179 ^He atoms, interacting through 
an accurate Aziz pair potential ' ^ in an almost cubic simulation box matching the 
periodicity of an hep lattice made up of A^, = 180 sites at a density p = 0.0294 
A"^. As usual, we apply periodic boundary conditions to the simulation box to 
emulate the infinite dimensions of the bulk system. 

Thanks to the accuracy of the Chin approximation for the action, it is possible 
to reach convergence of the physical observables in the limit e — !> with a rather 
large value of the imaginary time step e, making thus feasible the simulation of 
the quantum system with a small number of beads, even at low temperature ^\ 
More precisely, the convergence of the one-body density matrix pi (r) is achieved 
with a time step e = 0.033 K" ^ 

PIMC results for pi{r) at different temperatures and at the density quoted 
above are shown in Fig.[T] We have plotted, in the same figure, the zero-temperature 
estimation of pi for the same system and for a perfect hep crystal, obtained with 
the Path Integral Ground State method^'. We notice that the highest temperature 
at which the system presents a non-zero condensate fraction no, indicated by a 
plateau in the large r behavior of pi (r), is Tq = 0.2K. Our estimation for no at this 
temperature is mq = (8.4 ± 0.8) x 10~^. At temperatures T > Tq, the plateau dis- 
appears and the one-body density matrix presents an exponential decay at large r, 
which becomes more pronounced as the temperature increases, up to T = 0.75 K. 
Above this temperature, the decay of pi (r) becomes independent of T and it is 
similar to the large r behavior of pi (r) in commensurate (perfect) crystals. 

A relevant feature of PIMC is the possibility of giving a qualitative micro- 
scopic description of atomic positions in "^He crystals by means of its ring-polymer 
representation; the spreading of beads of each polymer gives an account of the 
delocalization of the atoms due to their zero-point motion. In Fig. |2] we show 
snapshots of the typical configurations of the incommensurate crystal at different 
temperature, plotting two-dimensional projections of the PIMC polymers lying in 
a basal plane of the hep lattice. 

At the highest temperature that we have studied, that is T = 1 K, the polymers 
do not spread on distances larger than the interatomic distance, indicating that the 
^He atoms tend to stay localized around their equilibrium positions. In this case, 
the vacancies are easily detectable inside the lattice. This behavior explains the 
fact that, at this temperature, the presence of vacancies does not affect noticeably 
the overall behavior of pi which, for the incommensurate crystal, is similar to the 
one of the perfect crystal. 

At lower temperatures, the delocalization of ^He atoms increases, allowing 
polymers to occupy different lattice points. In typical configurations at T = 0.5 K, 
we can detect polymers which spread on distances higher than the interatomic 
distance. This eventuality, which indicates that the vacancies begin to delocalize, 
enhances the possibility for the different polymers to superpose and thus to per- 
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Fig. 1 The one-body density matrix pi(r) for an hep crystal with vacancy concentration X,. = 
1/180 at density p = 0.0294 A^^^ and at different temperatures: T = IK (triangles down), T = 
0.75 K (triangles up), T = 0.5 K (diamonds), T = 0.25 K (squares) and T = 0.2K (circles). The 
dotted and dashed lines represent pi [r] at zero temperature respectively for the commensurate 
(X,, = 0) and incommensurate crystal (X,, = 1 / 1 80) at the same density, taken from Ref - ' . 



mute. However, at this temperature the polymers spreading on different lattice 
points are still rare events and the formation of long permutation cycles, which are 
necessary for BEC, is inhibited. 



If we further decrease the temperature, the polymers spreading on different 
lattice points become more frequent and it is impossible to associate a well de- 
fined lattice point to every quantum particle, as in the case at higher temperature. 
At r = 0.2K, the typical configuration of the incommensurate crystal looks like a 
commensurate system, indicating that the vacancies at this temperature are com- 
pletely delocalized and thus undetectable inside the crystal. The frequent occur- 
rence of a same lattice point occupied by beads belonging to different polymers 
strongly enhances the exchanges between the ^He atoms and allows also the cre- 
ation of long permutation cycles. In particular, for T < 0.2 K, it is possible to 
sample configurations presenting a non zero winding number, as the one shown in 
the low left panel of Fig.|2] In this picture, we show two following basal planes of 
the incommensurate hep lattice and we highlight a path which winds around the 
boundary conditions of the simulation box. The appearance of non-zero winding 
number paths in the sampled configurations is a clear signal that the simulated 
incommensurate crystal supports superfluidity at temperatures below To = 0.2K. 
On the other hand, this onset temperature is observed to be an increasing function 
of the vacancy concentration^^. 




Fig. 2 Two-dimensional projection of basal planes of the incommensurate hep crystal at differ- 
ent temperatures, represented according to the PIMC isomorphism of the classical polymers. At 
r = 1 K (high-left panel) the vacancy is localized and indicated by the circle. At T = 0.5 K 
(high-right panel), the vacancy begins to delocalize: the ellipse indicate a quantum particle delo- 
calized over two different lattice sites. AtT = 0.2 K (low-left panel), the vacancy is completely 
delocalized and cannot be easily detected. Below T = 0.2 K (low-right panel), the delocaliza- 
tion of vacancy strongly enhance the exchange between the bosons and allows the appearance 
of paths presenting a non-zero winding number, like the one represented by the thick line. 



4 Vacancies and the elastic constants 



We have carried out a computational study of the elastic properties of perfect (e.g. 
free of defects) and incommensurate solid ^He in the hep structure based on the 
diffusion Monte Carlo approach. This zero-temperature study is intended to im- 
prove our understanding of the response of solid helium to extemal strain, and 
extends the work initiated by Pessoa et al.^^ . In particular, we provide the zero- 
temperature dependence of C44 on pressure up to ^ 1 10 bar. This is a significantly 
higher pressure than previously considered both experimentally and theoretically. 
Our results are compared to experimental data and other calculations when avail- 
able and, as it will be shown later on, good agreement is generally found. The 
computational method that we employ is fully quantum and virtually exact, that 
is, in principle only affected by statistical uncertainties. In this sense, our study 
also represents an improvement with respect to previous zero-temperature first- 
principles work^^ based on variational Monte Carlo calculations (i.e., subject to 
bias stemming from the choice of the trial wavefunction). 



For small strains, the zero-temperature energy of a crystal can be expressed as 

1 '' 
E = Eq + -Vq ^ CijSiSj , (3) 

where Vq and Eq are the volume and intemal energy of the undistorted solid, {Qj} 
the elastic constants and {s,} the strain components defined such that ^i, ^2 and 
sj are fractional increases in the x, y and z directed axes, and S4, S5 and S(, angu- 
lar increases of the xy, xz and yz angles. ^^'^^ The symmetry of the crystal under 
consideration defines the number of elastic constants which are non-zero. In hep 
crystals, this number reduces to five namely Cn, C12, C33, C13 and C44, where C44 
is commonly known as the shear modulus. To calculate these elastic constants, is 
necessary to compute the second derivative of the intemal energy of the crystal 
with respect to the strain tensor a,y. For this, we express the primitive hep unit cell 
in terms of the translational vectors 



(4) 

(5) 

(6) 

where a and c are the lattice parameters in the basal plane and along the z axis 
respectively, and i, j and k correspond to the usual unitary Cartesian vectors, and 
two-atom basis set ri = ^ai + ^82 + fas and r2 = (0,0,0). 

In particular, the shear modulus C44 quantifies the response of the hep crystal 
to a heterogeneous strain in which the angle between by the c-axis and basal plane 
is tilted and the volume of the unit cell kept fixed. Such a shear deformation can 
be expressed as a transformation between sets of primitive translational vectors, 
{ai,a2,a3} -> {61,62,63} , where^ 





63 = ck , (7) 

e being a dimensionless parameter. It follows that 

^44 = 77- 3-7 , (8) 

where the equilibrium condition is fulfilled at e = 0. 

The simulation box used in our perfect (defected) pure shear calculations con- 
tains 200 (199) ^He atoms and was generated by replicating the hep unit cell 5 
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Fig. 3 Pressure-dependence of the calculated shear modulus of perfect and defected solid '^He 
at zero temperature. Experimental data are taken from Refs. ^* (Crepeau et al.),-- (Greywall), 
and^** (Syshchenko et al.) Linear tits to our results are also shown. 



times along ei and e2, and 4 times along 63. In proceeding so, hexagonal symme- 
try in our supercell calculations is guaranteed by construction. Periodic boundary 
conditions were imposed along the three directions defined by the edges of the 
non-orthorhombic simulation box. Prior to our shear modulus calculations, we 
determined the value of the equilibrium c/a ratio at each volume. We found that 
regardless of the pressure considered the optimal c/a value was always 1.62(1). 
In order to express our C44 (V) results as a function of pressure we employed the 
equations of state P{V) reported in Ref. ^^ (in both perfect and defected cases), 
which were deduced employing the DMC method and considering accurate finite- 
size corrections to the total energy. ^^ Additional details of our elastic constant 
calculations can be found in Ref. '^. 

In Fig. [3] we plot the pressure dependence of the shear modulus of perfect 
and defected solid ^He as obtained in our T = calculations. We found that 
these results can be accurately reproduced with a linear function of the form 
C44 (P) = 044 + b44P. In the perfect crystal, parameters 044 and Z?44 adopt the val- 
ues 92.2 (1.7) bar and 1.83 (0.04) respectively, whereas in the defected structure 
these are 93.3 (1.7) bar and 1.72 (0.04). It is observed that the elastic properties 
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C12 



C|3 



C33 



560(6) 

574(6) 
405(4) 



C44 



Commensurate 

Incommensurate 

Expt. 



210(3) 

227(3) 
213(4) 



130(2) 
152(2) 
105(13) 



639(7) 
649(2) 
554(22) 



140(2) 
138(2) 
124(2) 



Table 1 Elastic constants of the perfect (commensurate) and defected (incommensurate) hep 
solid '*He calculated with the DMC method at P = 26 bar Experimental values obtained by 
Crepeau et al— at a pressure P = 25.33 bar. Numbers within parenthesis stand for errors. 
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of perfect and defected helium are very similar at pressures close to melting (that 
is ^ 25 bar). Also we note that both sets of C44 data are compatible with exper- 
iments performed at low pressures and temperatures. ■^^'^LZ^ A similar agreement 
is achieved with the other elastic constants that we report at a pressure close to 
melting in Table 1 . At this particular pressure, all the constants of the defected 
crystal except C44 are larger than the ones of the perfect solid, with differences 
of variable size depending on the particular C,y. As compression is increased, the 
shear modulus of the defected crystal becomes appreciably smaller than that of 
perfect "^He. In relation to recent shear modulus experiments ^^ we must conclude 
that, based on this preliminary computational study, the increase of C44 observed 
at low-T can not be explained in terms of point defects. 



5 Summary and conclusions 

In the present work, we have reported recent results on the microscopic proper- 
ties of hep ^He with a tiny fraction of point defects, i.e., vacancies. To make our 
study as free of approximations as possible we have used state-of-the-art quan- 
tum Monte Carlo methods, both at zero and finite temperature. Relying only on 
the Hamiltonian, we have studied two particular aspects that we think relevant in 
the present discussion about the possible supersolid ^He phase. First, we have ad- 
dressed the question of the behavior of vacancies when the temperature decreases, 
approaching the temperatures where NCRI is observed. Our PIMC results unam- 
biguously show that there is an onset temperature for Bose-Einstein condensation 
that depends on the vacancy concentration. For the lowest concentrations that we 
can access with our approach this onset temperature is close to the experimental 
one. At this same temperature, the vacancies are observed to change their spatial 
behavior. Our simulations show that for T larger than the onset temperature, va- 
cancies appear as classical entities in the sense that we can always identify where 
they are, as it happens in a classical crystal. For temperatures smaller than the 
onset one, the vacancy is completely delocalized and the polymers representing 
the particles in the PIMC formalism spread over all the sites as if the crystal were 
perfect. 

In a second part, we have studied the elastic behavior of solid ^He in the limit 
of zero temperature. Technically, this is achieved by calculating energies of the 
deformed crystal by means of the DMC method. This approach is quite standard in 
classical simulations but much less worked out in quantum Monte Carlo. We have 
reported results for all the elastic constants at melting and the pressure dependence 
of C44. The presence of vacancies is also studied. Our preliminary results show a 
globally small effect with the more relevant result being the different slope of the 
elastic constants with the pressure: the slope for the incommensurate crystal is 
systematically smaller than the one for the commensurate phase. 

The role of vacancies in a supersolid scenario for solid ^He has been frequently 
excluded due to two arguments: the relatively high energy cost of its formation 
(above 10 K) and the possibility that vacancies form aggregates and eventually 
evaporate. It is certainly difficult that vacancies appear at very low temperature by 
thermal activation but we think that we cannot exclude the appearance of point 
defects along the growing process. Concerning the possible aggregation of vacan- 
cies due to a short-range attraction between them, there are different theoretical 
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predictions with contradictory observations. We are currently exploring this issue 
as a function of temperature. 
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